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The mobility of an overdamped particle, in a periodic potential tilted 
by a constant external field and moving in a medium with periodic friction 
coefficient is examined. When the potential and the friction coefficient have 
the same periodicity but have a phase difference, the mobility shows many in- 
teresting features as a function of the applied force, the temperature, etc. The 
mobility shows stochastic resonance even for constant applied force, an issue of 
much recent interest. The mobility also exhibits a resonance like phenomenon 
as a function of the field strength and noise induced slowing down of the par- 
ticle in an appropriate parameter regime. 
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I. INTRODUCTION 

The simple phenomenon of stochastic resonance ( SR) has been an active field of research 
for more than one and a half decades. Its appeal to natural processes has generated a lot of interest 
among workers from various branches of science and engineering JTJ . This interest is largely because 
of the fact that in this phenomenon noise plays a useful role. The presence of noise helps in improving 
the quality of response of a system subject to an alternating external field. The conventional SR is all 
about the( output signal-to- noise ratio) optimization of an input signal by noisy nonlinear systems 
by suitably tuning the noise level. It is reflected as a peak in the output signal-to-noise ratio as a 
function of input noise strength. To observe SR one requires three basic ingredients (i) an energetic 
activation barrier, (ii) a weak coherent input periodic signal, and (iii) a source of noise. Attempts 
are being made to reduce the number of constraints for the realization of SR. In our present work, 
we are interested in the response optimization, in the presence of a non-oscillating external signal. 

There was a suggestion that SR could be observed in the drift velocity of an overdamped 
Brownian particle in a tilted periodic potential 0. The suggestion was an important one for it 
sought to do away with the oscillating input signal, considered to be an essential ingredient for 
the observability of the phenomenon of SR, and replace it with a constant uniform external field. 
Though the suggestion has now been proved to be incorrect |||4]] for overdamped particles, it is 
argued to be correct for underdamped particles || where the inertia of the particles could act as 
surrogate to the external oscillating field in a periodic potential. The present work shows that SR 
can indeed be observed in the drift velocity of overdamped particles in a periodic potential subjected 
to non-oscillating uniform constant field but in a medium where the particles experience spatially 
periodic frictional drag. 

The variation of space dependent friction 77(g) ( in the absence of spatial variations in 
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temperature) influences the dynamics of the particles in a potential field and helps the system to 
approach towards its equilibrium or steady state. The relative stability of the competing states 
is governed by the individual Boltzmann factor, i.e inhomogeneity in friction does not affect the 
stationary ( equilibrium) properties of a system but it does affect the dynamical ( nonequilibrium) 
properties such as the relaxation rates. In contrast, temperature inhomogeneity, changes the relative 
stability of the otherwise locally stable states ( thereby creating new steady states). Moreover, the 
relative stability of the system in two different local minima depends sensitively on the temperature 
profile along the entire pathway connecting the two local minima. Temperature nonuniformity 
together with spatially varying friction coefficient or alone has been shown to give rise to a net 
current in a periodic potential P^TT| in the absence of an applied bias. Spatially varying friction 
coefficient, 77(g) of the medium, alone, however cannot give the net probability current j in such a 
potential. But when, in addition, a constant force F is applied a finite net current may result at a 
constant temperature T of the medium. We consider a sinusoidal potential with period 2tt and a 
sinusoidally varying model friction coefficient of the same period but with a phase difference <fr and 
modulating amplitude A. The mobility H|12| /j,(=<-4>/F= ^P, where < -4 > is the average 
drift velocity of the particles) shows very interesting behaviour in the parameter space of (T, F, A, 4>). 
Apart from the above mentioned SR as a function of noise, the mobility also exhibits resonance like 
phenomenon as a function of F and the phenomenon of noise-induced stability etc. 

Before we present the results of our calculation in section 3 we discuss the necessary for- 
malism in section 2. We summarize our work in the last section 4, with discussion. 

II. THE MOBILITY 

The motion of an overdamped particle, in a potential V(q) and subject to a space dependent 
friction coefficient 77(g) and an additional constant force field F, at temperature T, is described by 

3 



the Langevin equation P, |I0| , |I2| . |I"3 



dq 
~dt 



(V\q) - F) 



(1) 



77(g) 7](q) 2 V 77(g) 

Here e(t) is the randomly fluctuating Gaussian noise term with the property 

< e (t) >= , 
< e(t)e(ti) >=25(t-t l ) , 

where < .. > denotes average over the distribution of the fluctuating quantity e(t). The primes in 
Eqn. (1) denote derivative with respect to the space variable q. Eqn. (1) has been derived earlier 
using microscopic treatment of system-bath coupling |f3~0| , [T3[1 . It should be noted that the above 
equation involves a multiplicative noise with an additional temperature dependent drift term. The 
additional term turns out to be essential in order for the system to approach the correct thermal 
equilibrium state. The motion is equivalently described by the Fokker-Planck equation [fT|, p~Q| , [T3| , [T^ 



dP( q ,t)_d 1 .^a^fl + n 



(2) 



dt dq r](q) 1 a dq 

The stationary probability current j, for periodic V(q) and 77(g) with the same periodicity 2ir, is 
given as fl6l|T2| 

k B T(l - 



-2ttF 

7r' 



J 



„ -V(y)+Fy 



dy Jy * r](x)e k B T dx 



(3) 



and the corresponding mobility P,|T2fl jj, = Clearly j — > as F — > 0, but jj, remains finite for 



finite temperature T. The mobility in this linear response regime []T2| is given by 

hm /i = -yj-. ^~J. 



(4) 



Notice that the above expression for /1 involves a combination of rj(x) and V(x) which contribute in 
different ways and cannot be accounted for by a single "effective" potential. We would like to point 
out that a particle in a medium, with space dependent temperature profile, experiences an effective 



potential ( or net force). This in turn, is able to shift the stable points and introduces noise induced 
transitions ||. In our present work we emphasize that the space dependent friction influences the 
kinetic properties of the system and unlike in the case of space-dependent temperature we do not 
have a net additional force in the steady state. 

In the absence of any external periodic potential the drift velocity is linear in F and the 
corresponding mobility is given by /i — 4, where fj = ^ Jy +27T r](x)dx } is the average value of 
the friction coefficient over a period. The above result follows from straight forward mathematical 
analysis. From Eqn. (3) one can see that the magnitude of // depends sensitively on the potential and 
the frictional profile over the entire period. And it can be much larger or smaller than the mobility 
of a particle moving in a homogeneous medium with friction coefficient equal to fj. Depending on the 
temperature and the applied field, mobility can in fact be much larger than the asymptotic mobility 
( with respect to F and T). This implies that at some particular values of F and T the mobility, 
in the absence of periodic potential is much smaller than in the presence of periodic potential. This 
counter-intuitive result emerges because of the complex manner in which the phase shift between 
the potential and the frictional profile, the amplitude of the friction coefficient, etc, influence the 
mobility. We consider for our calculation, 

V(q) = — sin q, (5) 
■q(q) = 770(1 - Asin(g + </>)), (6) 

with < A < 1 so that fj = tjq. One can easily see that for F < 1, as T — > 0, \i — > 0. Also, as 
F, T — > 00, \i — > Moreover, j(F) ^ —j(—F), except when = and tt, and hence mobility is 
asymmetric with respect to the force. For intermediate values of T and F one needs to evaluate the 
double integral in the denominator of Eqn. (3) numerically ||15||and explore the variation of \x : 



We calculate fi in various sections of the parameter space of (T, F, A, 0). Henceforth we have scaled 
mobility, temperature and force to dimensionless units as in Ref. |[L2|| . 



III. RESULTS. 

The mobility fx shows many interesting features even at T = 0. At T = 0, following Ref. 
fl"2|| closely, fx remains zero for \F\ < 1 for all values of and A. However, and A play important 
role for \F\ > 1 and the mobility is given by 

in _ j_w2 



i + A(l-^/l-^)sin0 
For applied fields smaller than the critical field (F = 1) the particle remains in a locked state or 
confined to the local minima of the potential. For F > 1 (i.e, when the potential barrier for motion 
does not exist) the particle moves down the potential slope or will be found in a running state P,|12~ 



Observe that at T = 0, ix(F, 0) = fx(—F, —0) and this happens to be true for T ^ as well. We 
may roughly divide the full range of (0-27r) into two regions where fx shows qualitatively distinct 
features [Fig. (1)]. 

From Fig. (la) we observe that for F > 1, < < 7T, rj Q fx always remain less than one but 
approaches one asymptotically as F — > oo. Fig. (lb) shows that in the range 7r < < 27r, rjofi could 
even be larger than one for large A and exhibits a maximum. This occurrence of maxima is a new 
feature and could only be ascribed to the space dependence of friction coefficient T](q). It is to be 
noted that the value of the mobility at the maxima is always larger than the asymptotic mobility 
(4). In Fig. (2) we have given a surface plot of the mobility obtained from Eqn. (8) as a function of 
F and 0. This figure clearly shows that there is a small region in 0, around 1.27T where the maxima 
can be observed, elsewhere mobility monotonically increases with F. In contrast, in a medium with 
a space independent friction , the mobility in the high friction limit monotonically increases with F 



for all values of T and asymptotically goes to one [I2|. For given F > 1, however from Eqn. (8), 
maximum occurs at = -y- and minimum at <f) — |. It should be noted, that the space dependence 
of the friction coefficient does not alter the threshold value (F = 1) for nonzero mobility. 

As the temperature T is increased from zero, the thermal fluctuations make the current 
j nonzero even for \F\ < 1, when the barrier to free passage motion is nonzero. For T ^ 0, as 
F — > 0, j ' ^ linearly with F but the mobility remains finite [Eqn. (4)]. For a given temperature 
T 7^ and when the field is very high, the effect of the periodic potential vanishes and consequently 
rjQ^i — > 1. However in the small F regime 770/i can be much larger than one, leading to a peak as a 
function of F depending on the parameter values of 0, A,T. The resonance feature as a function of 
F can also be observed at finite temperatures, as shown in Fig. (3). 

Let us now examine the variation of 770/i as a function of the modulation parameter A of 
the friction coefficient. Fig. (4) shows the variation of 770/i for = at fixed T = 2.0 and Fig. (5) 
for = 7r at the same T ( For = and n, 770/i is symmetric about F = 0). When V(q) and 77(g) 
are in phase, the mobility monotonically decreases with A (F = having the least value at A = 1 ) 
whereas when they are in opposite phase mobility increases with A. For = tt and large A ( Fig. 
[5]), 770/i becomes larger than one for all values of F. Figures (4 and 5) essentially corroborate the 



observations made in Ref. |Uj that when = the friction coefficient is the largest at the positions 
where the particle has the largest acquired velocity due to the potential, and it is smallest where the 
particle has the smallest velocity. Thus the damping term has the severest effect of slowing down 
the particle. When = tc the situation is just the opposite. In this case the frictional effect is the 
least. In fact it has a positive effect and as can be seen from Fig. (5) the mobility 770/i even becomes 
larger than one. For any other 0, however, the situation is complicated to analyze as simply. To 
emphasize this we have given a surface plot of 770/i with A and in Fig. (6) for T = 2.0 and F = 0.8. 
Mobility increases monotonically with A in a region for around I.I71", and elsewhere it decreases 



monotonically. 

We now discuss the variation of rjo/i as a function of temperature T ( or the noise strength). 
At this point it is pertinent to note that there has been discussion in literature about the motion of 
an underdamped particle, in tilted ( due to a constant force F) but otherwise periodic potentials, 
being in two states depending on the value of F and the damping constant 770 of the medium ||[12 



The particle could be in a locked state ( remaining confined to a local minimum of the potential) 
or in a running state ( corresponding to a motion down the potential slope). The two states could 
also coexist even for \F\ < 1 ( for periodic potential V(q) as given in Eqn. (5) ). However, in the 
overdamped case at T = the particle can be either in the locked state (for \F\ < 1) or in the 
running state( for \F\ > 1). But as T is increased the particle incoherently hops from one stable 
( locked state) to another one even for \F\ < 1, leading to a net current. When in addition, we 
consider the friction coefficient to be space dependent, the mobility as a function of temperature 
shows many more interesting features depending on the values of F, A and 0. We confine our attention 
to large values of A ( but < 1), for example, A = 0.9. 

In the entire range of 0[O < < 2tt] we present results for a few typical values of 0. 
Fig. (7a) shows the nature of SR at = 0.97T for arbitrarily selected values of F. We see that for 
some values of F, especially for F < 1, 770// becomes larger than one at intermediate values of T 
( the peaks appear at larger T for larger F) and asymptotically tends to one as T — > 00, for all 
values of F. We find SR even for F — > for F — > and in fact the peak mobility is larger for 
F < 1 than for F > 1. Similar behaviour could be seen for < 0.9n also. These peaks for small 0, 
however, are very broad. Fig. (7b) is for = 1.447T. In this figure the peaks are sharper for larger 
values of F (F > 1) and the peaks disappear for F < 1. The temperature T(= Tp) corresponding 
to the peak value of mobility is much lower than that in the previous case (0 = 0.97r) and unlike the 
previous case it decreases and tends to zero as F is increased. Figs. (7a) and (7b) clearly exhibits 



the occurrence of SR for F < 1 as well as for F > 1. It should be noted that we observe stochastic 
resonance for smaller values of A(< 0.9) too. But as A is decreased the peaks become less prominent 
and the range of over which the phenomena could be observed shrinks and finally disappears. 

We now turn to examine the location of maxima Tp, as a function of system parameters. 
In Fig. (8) we have plotted Tp against for various values of F, the inset shows the variation of Tp 
versus F for some typical values of as indicated on the figure. For small fields F < 1, Tp versus 
is a non-monotonic function exhibiting a minima, which shift towards higher values of as F as 
increased. As we go away from the minima on both sides, the resonance behaviour becomes more 
broader. For beyond 1.57T, mobility does not show a maxima( or Tp). For values of field larger than 
about 1.2, Tp decreases monotonically with 0. Tp shows a monotonic or a non-monotonic behaviour 
as a function of F as can be seen from the inset of Fig. (8). This behaviour is very sensitive to 
material parameters. In Fig. (9) we have given a surface plot of the mobility as a function of 
and T for F = 0.5 and A = 0.9. For values ranging from = 0.97T to 1.44tt we obtain SR in the 
mobility as a function of temperature. It is to be noted that this surface plot as well as all the earlier 
surface plots exhibit a maxima in the mobility as a function of for a given temperature T and field 
strength F, on the other hand the mobility shows a monotonic behaviour with A (see for example 
Fig. (6)). This seems to suggest that may play the characteristic role of frequency in our model 
in the absence of additional external frequency signal P,[TT| . 



Fig. (10a) is for = l.bir. For forces less than a certain F (closer to the critical value of 
F — 1) all curves begin with rjofi less than one and gradually increase with T and asymptotically 
reach the value one. And for larger forces the mobility always remain larger than rjofi = 1 and 
asymptotically reach one as T —>■ oo. For large values of F(> 1), for which the potential barrier for 
motion has already disappeared, the mobility decreases as T is increased from T = 0. This implies 
that the motion of the overdamped particle becomes more sluggish than at T = 0, where one would 



have expected the particle to become more mobile as T is increased from zero. Our observation 
indicates that the presence of noise actually slows down the motion of deterministically overall 
unstable states in the appropriate range of F and cf). This is somewhat akin to the phenomenon 



of noise-induced stability of unstable states 18-21]. For the values of F larger than the critical 



field, the barrier to the motion of the particle disappears and in this case ( depending on the system 
parameters) the noise slows down the motion between unstable states whereas for F < 1 the noise 
helps to overcome the potential barrier for the passage of the particle. And finally in Fig. (10b), 
for = 1.67T , we observe some very interesting features that are somewhat contrary to common 
expectations. Values of F very close to ( but larger than the critical value one), the mobility attains 
a minimum value and then begins to increase with T and finally reaches tjq/j, = 1, asymptotically, as 
T — > oo. In this case we observe a cross-over from the regime of noise-induced slowing down to the 
noise-enhanced mobility for a given F as a function of the temperature (noise strength). All these 
interesting features ( some of them even counterintuitive) described above, results from a subtle 
combined effect of the periodic space dependent friction and the periodic potential in the presence 
of a constant applied force. The phase difference <fi between the potential and the friction coefficient 
plays a crucial role in determining their nature. 

IV. SUMMARY AND DISCUSSION 

We have theoretically studied the motion of an overdamped particle in a tilted( by a constant 
force) but otherwise sinusoidal potential in an inhomogeneous medium. We considered the particle 
to experience a sinusoidally varying friction coefficient while in motion, but the temperature of the 
medium is kept uniform . The behaviour of the motion of the particle characterized by its mobility, 
is influenced by the phase difference between the potential function and the friction coefficient in 
important ways. If the friction coefficient is small where the potential is high, then the bias F will 



in 



influence the mobility differently than when the friction coefficient is large where potential is high. 
By choosing the phase difference suitably the mobility could be enhanced by tuning the noise level 
( temperature) of the system or lowered depending on the value of the phase difference. The former 
indicating the occurrence of stochastic resonance in the motion of the overdamped particle in a 
periodic potential with the application of a non-oscillating constant field. The decrease of mobility, 
on the other hand, indicates that the motion could be made sluggish by enhancing the noise level 
( temperature). We have also observed the peaking behaviour in the mobility as a function of applied 
field. Of course, both these effects cannot be seen in a homogeneous medium with constant friction 
coefficient and uniform temperature. 

At this point it is worthwhile to note that in the conventional treatment on SR it has been 
established that SR occurs when Kramer's time of barrier passage roughly matches the ( dominant 
characteristics) time scale of the input ( externally applied) signal. But in the present work there is 
no time varying externally applied input signal and hence the criterion for the occurrence of SR is 
not satisfied. However some recent works have shown that the SR can also occur in the absence of 
an external periodic force, which is a consequence of the intrinsic dynamics of the nonlinear system 



22| , |23|1 . Such a behaviour has been described as autonomous SR. It may also be possible that the 
SR phenomena that we have discussed, may be a little closer to dither( or threshold SR) . 

In this system, that we study, the current is in general asymmetric with respect to the 
reversal of the direction of the applied force. This result has earlier been used to obtain a net 
unidirectional motion in the presence of time oscillating field, but with net average force zero( i.e 
without an obvious bias) 
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FIGURE CAPTIONS 



FIG. 1. Mobility i] fi as a function of F at T = for various values of A (all in dimensionless 
units ) for (a) = 0.9tt, (b) = 1.2tt. 

FIG. 2. Surface plot of Mobility at T = as a function of F and at A = 0.9. is in units 

of 2tt. 

Fig. 3. Surface plot of mobility rjopL as a function of applied bias F and at T = 0.8 and 
A = 0.9. is in units of 2-n. 

FIG. 4. Mobility i] fi as a function of A at T = 2.0 and = for various values of F . 

FIG. 5. Mobility i] fi as a function of A at T = 2.0 and = tx for various values of F . 

FIG. 6. Surface mobility as a function of A and at F = 0.5 and T = 2.0, to highlight the 
monotonic behaviour of mobility with A. 

FIG. 7. Mobility r] fi as a function of temperature T at A = 0.9 for various values of F for 
(a) = 0.97T (b) = 1.447T . The insets highlight the maxima in the curves. 

Fig. 8. Temperature corresponding to the peak mobility (T P ) as a function of ( in units 
of 2n ) for various values of F. The inset shows the variation of Tp as a function of F for various 
values of 0. 

Fig. 9. Surface plot of mobility rjo/i as a function of temperature T and at F = 0.5 and 
A = 0.9. is in units of 2n. 

FIG. 10. Mobility rjo/j, as a function of temperature T at A = 0.9 for various values of F for 
(a) = 1.57T (b)0 = 1.67T. The inset of figure (b) highlights the minima in mobility. 
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